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We investigate adding two sterile neutrinos to resolve the apparent tension existing between short- 
baseline neutrino oscillation results and CPT- conserving, four-neutrino oscillation models. For both 
(3+1) and (3+2) models, the level of statistical compatibility between the combined dataset from 
the null short-baseline experiments Bugey, CHOOZ, CCFR84, CDHS, KARMEN, and NOMAD, 
on the one hand; and the LSND dataset, on the other, is computed. A combined analysis of all 
seven short-baseline experiments, including LSND, is also performed, to obtain the favored regions 
in neutrino mass and mixing parameter space for both models. Finally, four statistical tests to 
compare the (3+1) and the (3+2) hypotheses are discussed. All tests show that (3+2) models fit 
the existing short-baseline data significantly better than (3+1) models. 
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There currently exist three experimental signatures for 
neutrino oscillations. The two signatures seen originally 
in solar and atmospheric neutrinos have been verified by 
several experiments, including experiments carried out 
with accelerator and nuclear reactor sources. The results 
on atmospheric neutrinos can be explained by I'u dis- 



appearance due to oscillations Q, 0, |3j , while those on 
solar neutrinos can be explained by oscillations 0, . 
The third signature is Df, appearance in a beam, ob- 
served by the short-baseline, accelerator-based LSND ex- 
periment at Los Alamos 'oj . This signature is strong from 
a statistical point of view, being a 3.8(t excess, but further 
confirmation by an independent experiment is necessary. 
The MiniBooNE experiment at Fermilab will be able to 
confirm or refute the LSND signature in the near future, 
with an experimental setup providing different systemat- 
ics and higher statistics than LSND Q. 

Taken at face value, the three experimental signatures 
point to three independent mass splittings. Three neu- 
trino masses do not appear to be able to explain all of 
the three signatures (see, however, 0|). One way 

to solve this puzzle is to introduce different mass spectra 
for the neutrino and antincutrino sector, thereby requir- 
ing CPT- violation but no extra neutrino generations 
Another possibility is to add additional neutrinos with no 
standard weak couplings, often called "sterile neutrinos" . 

In this paper we assume CPT- and CP-invariance, and 
we explore the possibility of adding one or two neutrino 
generations beyond the three active flavors assumed by 
the Standard Model. We focus on extensions of the neu- 
trino sector where the addition of fourth and fifth mass 
eigenstates are responsible for the high Am? LSND oscil- 
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FIG. 1: Flavor content of neutrino mass eigenstates in (3+1) 
models (a), and (3+2) models (b). Neutrino masses in- 
crease from bottom to top. The Ve fractions are indicated by 
right-leaning hatches, the Vfj, fractions by left-leaning hatches, 
the u-r fractions by crosshatches, and the Vs fractions by no 
hatches. The fiavor contents shown are schematic only. 



lations, and the three lower mass states explain solar and 
atmospheric oscillations. When only one sterile neutrino 
is added, these models are labelled as (3+1). The fiavor 
content of the four neutrino mass eigenstates for these 
models is schematically shown in Fig. The (3+1) hi- 
erarchy in Fig. nji is as opposed to the (2+2) hierarchy, 
where the solar and atmospheric mass splittings are sep- 
arated from each other by the LSND Am^. The (2+2) 
models require a different global analysis from the one 
discussed in this paper. The simplest (2+2) models ap- 
pear to be only marg inally consistent with neutrino os- 
cillations data H, [ia |. even though more general (2+2) 
mass and mixing scenarios might represent a viable solu- 
tion to explain solar, atmospheric, and LSND oscillations 

The (3+1) models are motivated by the criterion 
of simplicity in physics, introducing the most minimal 
extension to the Standard Model that explains the ex- 
perimental evidence. However, theories invoking ster- 
ile neutrinos to explain the origin of neutrino masses do 
not necessarily require only one sterile neutrino. Indeed, 
many popular realizations of the see-saw mechanism in- 
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troduce three right-handed neutrino fields [Illlllll. In 
particular, (3+2) neutrino mass and mixing models can 
be obtained in several see-saw mechanisms From 
the phenomenological point of view, it is our opinion 
that two- and three-stcrilc neutrino models should also 
be considered and confronted with existing experimen- 
tal results. In this paper, we consider the results from 
the short-baseline experiments Bugev ITsIl , CCFR84 , 
CDHS IliiCHOOZ 113, KARMEN 22], LSND Q, and 
NOMAD|23l|. and examine how well (3+1) and (3+2) 
models agree with data. A schematic diagram for (3+2) 
models is shown in Fig. We do not consider (3+3) 
models in this paper. From our initial studies, we believe 
that the phenomenology of a (3+3) model is similar to a 
(3+2) model. 

The paper is organized as follows. In Section m we 
specify the neutrino oscillations formalism used in this 
analysis to describe (3+1) and (3+2), short-baseline, os- 
cillations. In Section ITTll and Hvl we present the results 
obtained for the (3+1) and (3+2) models, respectively. 
For both models, we first derive the level of compatibil- 
ity between the null short-baseline (NSBL) experiments 
and LSND. Second, we perform a combined analysis of 
all seven short-baseline experiments (including LSND) to 
derive the preferred regions in neutrino mass and mixing 
parameter space. In Section V, we discuss four statisti- 
cal tests to compare the (3+1) and (3+2) hypotheses. In 
Section VI, we briefly mention other experimental con- 
straints on (3+1) and (3+2) models. In Appendix 1X1 we 
describe the physics and statistical assumptions used in 
the analysis to describe the short-baseline experiments. 

II. NEUTRINO OSCILLATIONS FORMALISM 

Under the assumptions of CP- and CPT-invariance, 
the probability for a neutrino, produced with flavor a 
and energy E, to be detected as a neutrino of flavor /3 
after traveUing a distance L, is ^2^]: 

n 
j>i 

where a — e,^,T,s (s being the sterile flavor); U is the 
unitary neutrino mixing matrix; Xji = 1.27 Arrij^L/E; 
Am?i = m3 — mf', and n is the number of neutrino 
generations. Neglecting CP-violating phases, there are 
in general (n — 1) independent mass splittings, and 
— n — n{n ~ l)/2 independent mixing matrix ele- 
ments. The situation simplifies considerably by consid- 
ering short-baseline (SBL) data only. In this case, it is a 



good approximation to assume X21 = X32 = 0, and only 
(n — 3) independent mass splittings are present. More- 
over, given the set of SBL experiments considered, the 
number of mixing matrix elements probed is only 2(n— 3), 
as we show now for the (3+1) and (3+2) cases. 

For (3+1) models, n—A, and only one mass splitting 
Am^ = Amli ~ Am|2 — Aml^ appears in the oscil- 
lation formula: this is sometimes referred to as to the 
"quasi-two-neutrino approximation" , or "one mass scale 
dominance" [l^. Using the unitarity properties of the 
mixing matrix, we can rewrite Eq.^for (3+1) models in 
a more convenient way: 

(2) 

which depends on the mass splitting (Amli) and mix- 
ing parameters {Uai, U/34) of the fourth generation only. 
Since the two-neutrino approximation is satisfied in the 
(3+1) case, we can express Eq. |5|in the usual forms: 

P{va — > i^/3) = sin^ 26ai3 sin^ X41, 01 ^ (3 (3) 



P{l'a Va) = 1 ~ Sin^ 29aa sin^ Xii (4) 

where Eq. |31 applies to an oscillation appearance mea- 
surement, Eq. ^to a disappearance measurement. 

In this paper, we use the data from the Bugey, 
CCFR84, CDHS, CHOOZ, KARMEN, LSND, and NO- 
MAD experiments. Bugey and CHOOZ data constrain 
Ve disappearance; CCFR84 and CDHS data constrain 
disappearance; and KARMEN, LSND, and NOMAD 
data constrain — > v,, oscillations. Therefore, from 
Eas. I2I3I and^I the experiments constrain the following 
combinations of (3+1) mixing parameters: 

• Bugey, CHOOZ: sin^ 26*^6 = 4C/2^(1 - C/^J; 

• CCFR84, CDHS: sin^ 261^^ = AUl^{l - U^^^); 

• KARMEN, LSND, NOMAD: sin^ 26'^e = 'iU^4U^4- 

In (3+1) models, the tension between the experimental 
results comes about because Bugey, CHOOZ, CCFR84, 
CDHS, KARMEN, and NOMAD hmit the two indepen- 
dent mixing matrix parameters Uei and to be small, 
whereas LSND demands nonzero values. 

In (3+2) models, we introduce two sterile neutrinos. 
Using Eq. ^ ^-nd the unitarity of the mixing matrix, 
the (3+2) neutrino oscillation probability formula can 
be written: 



P(i/q Vp) = Saf} - A[{5al3 - UaAUpi - U abU f3^){U aiU pAS.i'C? X41 + UabUpz S\T? X51) + U aiU abU PaU p^ SIT? X54] = 

= <5q/3 + 4[;7^4([/|4 - Sap) Xii + Ul^iUj^ - Sap) siu^ X51 + UaiUpiUabUpzisin^ X41 + sin^ X51 - sin^ 2:54)] (5) 
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which in our case depends on two independent mass 
sphttings (Am|j^, Am^j^) and four independent mixing 
matrix parameters {Uai, Uab, with a = e, fi). Eq.|21can 
be recovered from Eq. |31 by requiring Ua5 = = 0. 
In (3+2) models, the quasi-two-neutrino-approximation 
is not vahd, since there are three distinct Am? values 
contributing in the oscillation formula: Am|j^, Arrigj^, 
and Am54, and therefore three distinct oscillation am- 
plitudes: (sin^ 26'q,/5)4i, (sin^ 20ct/3)5i, and {sin^ 29a|3)54■ 
We now comment on the Monte Carlo method 
used to apply the above oscillation formalism to 
the analyses presented in this paper. We require 
the neutrino mass splittings to be in the range 
0.1 eV^ < Amli, Am^i < 100 eV^ with Am^i > Amli. 
Each mass splitting range is analyzed over a 200 point 
grid, uniformly in logj^gAm^. In Section IVII we briefly 
discuss why large mass splittings are not necessarily in 
contradiction with cosmological (and other) data. The 
values of the mixing parameters, J7e4, U^i, C/es, and 
are randomly generated over a four-dimensional 
space satisfying the four requirements: U^^ + U^^ < 0.5, 
+ U^5 < 0-5: where: i = a ~ e,fi. These 

four inequalities are introduced to account for the fact 
that large electron and muon flavor fractions in the 
fourth and fifth mass eigenstates are not allowed by 
solar and atmospheric neutrino data. In principle, since 
the CDHS constraint on disappearance vanishes for 
Am^ ~ 0.3 eV^, as shown in Appendix 1X1 the upper 
limit on disappearance from atmospheric neutrino 
experiments above the atmospheric Arri^ should be 
considered instead. In this paper, we do not reconstruct 
the likelihood for atmospheric data that would give 
the exclusion region for i/^ disappearance in the range 
Am^t„j < Am^ < 0.3 eV^. However, the effect that 
the atmospheric constraints would have on our results is 
expected to be small. For example, in Ref. Bilenky 
et al. use the atmospheric up-down asymmetry to 
derive the upper hmit U^^^ < 0.55 at 90% CL for (3+1) 
models, which is satisfied by our initial requirements 
Ue4 + UnA < 0.5, + < 0.5. A more recent 
analysis |34| of atmospheric neutrino data using the full 
zenith angle distribution provides a tighter constraint 
on sin^ than the one given in Ref. |3^ ; the impact 
of this additional constraint on our SBL analysis is 
discussed in Sections Hill and IVll Finally, from Eqs.[21 it 
is clear that the relative sign of both Uei and t/p4 cannot 
be inferred in (3+1) oscillations. Similarly, from Eq. [S] 
the only physically observable relative sign between 
mixing parameters in CP-conserving (3+2) models is 
sign{Ue4U^4UebU^^)\ therefore, this is the only sign 
related to mixing parameters that we randomly generate 
in the analysis. 

Throughout the paper, we make use of the Gaussian 
approximation in determining allowed regions in pa- 
rameter space. In general, this means that the regions 
of quoted confidence level are the ones enclosed by 
contours of constant values, whose differences with 
respect to the best-fit value depend on the number 



of free parameters in the model 31]. In the text, we 
use the symbol 5 to denote the values of the confidence 
levels derived in this way. As pointed out in [3^ . 
this approach should be considered approximate, as it 
may provide regions in parameter space of both higher 
and lower confidence than the one quoted. Regions of 
higher confidence than the quoted value may result from 
the presence of highly correlated parameters. Regions 
of lower confidence may result from the presence of 
fast oscillatory behavior of the oscillation probability 
formula, Eq. 



III. RESULTS FOR (3+1) MODELS 

This section, like the next one on (3+2) models, con- 
sists of two parts. First, we quantify the statistical com- 
patibility between the NSBL and LSND results, follow- 
ing a method described in |2^ 0| , originally proposed to 
establish the compatibility between the LSND and KAR- 
MEN results. Second, we perform a combined analysis 
of the NSBL and LSND datasets, to obtain the favored 
regions in neutrino mass and mixing parameter space. 



A. Statistical compatibility between NSBL and 
LSND 

Many analyses of the NSBL experiments within (3+1) 
models have concluded that the allowed LSND region is 
largely excluded |2^|2^|33- Here, we repeat this study 
with two purposes. First, we use this study to give con- 
text to our discussion of the basic model and techniques 
which will be expanded in later sections. Second, we 
demonstrate that our fit, which forms the basis of our new 
results for (3+2) models, reproduces the expected (3+1) 
exclusion region. For a discussion of the physics and sta- 
tistical assumptions used to describe the short-baseline 
experiments used in the analysis, the reader should refer 
to Appendix 1X1 

In this section, the NSBL and LSND datasets are 
analyzed separately, providing two independent allowed 
regions in (sin^ 20^e, Am^) space. The level of statisti- 
cal compatibility between the two datasets can be deter- 
mined by studying to what degree the two allowed regions 
overlap, as will be quantified later in this section. 

For each randomly generated (3+1) model, we calcu- 
late the values for the functions X^sbl ^'^^ Xlsndi 
where x^ssi, is defined as: 

XnSBL — XBugcy^~XcHOOZ^~XcCFR84^~XcDHS^"XKARMEN^~XN0^4AD 

(6) 

For the analysis described in this section, the NSBL 
and LSND allowed regions are obtained using two differ- 
ent algorithms, reflecting the fact that the NSBL dataset 
provides upper limits on oscillations, while the LSND 
dataset points to non-null oscillations. 
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FIG. 2: Compatibility between the NSBL and LSND datasets in (3+1) models. Fig.|5li shows the 90% (grey solid line), 95% 
(black dotted line), and 99% (black solid line) CL exclusion curves in (sin^ 26l^e, Am^) space for (3+1) models, considering 
the null short-baseline (NSBL) experiments Bugey, CCFR84, CDHS, CHOOZ, KARMEN, and NOMAD. Fig. Ill also shows 
the 90%, 95%, and 99% CL allowed regions by our analysis of LSND data. Fig. Ha) is as Fig. H, but in (plsnd, Am^) space, 
where Plsnd is the LSND oscillation probability (see text for the definition). Fig. H) and d) show the minimum values as 
a function of Am^ for the NSBL and LSND datasets (143 and 3 d.o.f., respectively). 



The NSBL allowed regions at various confidence lev- 
els i5nsbl are obtained via a raster scan algorithm |32|. 
Let Xnsbl be the value for the particular model and 
(XNSBL)min Am2 be the minimum for the Arn^ value 
considered. For example, our quoted 95% CL upper limit 
on sin^ 20^e is given by the maximum value for the prod- 
uct 4[/g4[/^4 chosen among the models vifhich satisfy the 
inequality xLbl - (XNSBL)min,Am2 < 5.99. The value of 
5.99 units of x^ is chosen because there are two free pa- 
rameters Ue4, Ufj_4 for (3+1) models with fixed Am^. We 
note that even for the NSBL dataset, the parameters Ue4, 
U^A can be correlated, since the KARMEN and NOMAD 
results probe a combination of the two parameters. 

The LSND allowed regions at various confidence lev- 
els ^LSND are obtained via a global scan algorithm [s^ . 
For example, for ^lsnd — 0.95 we require xLnd ~ 
(XLSND)min < 5.99, where (xLSND)min is now the global 
LSND x^ minimum value, considering all possible Am? 
values. The LSND allowed region is computed for two 
free parameters as for the NSBL case, but the parame- 
ters are now Am? and U^^Uei, as opposed to U^^ and 
f/e4- Compared to the NSBL case, the number of free pa- 
rameters is reduced by one because the LSND D^^ De 
search only probes the product Uf^^Uei and not the two 
mixing matrix elements individually, and it is increased 
by one because the allowed region is now obtained by 
scanning over all possible Am^ values. 

The regions allowed in (sin^ 20pe, Am^) parameter 
space by both the NSBL and LSND datasets are shown 
in Fig. Hl. The NSBL allowed regions shown are 2- 



dimensional projections of 3-dimensional allowed regions 
in (Am^, Ue4, Ufi4) space. The NSBL results alone allow 
the regions to the left of the solid grey, dotted black, 
and solid black lines in the Fig. at a confidence level 
5^33^ = 0.90, 0.95, 0.99, respectively. In Fig. the 
Slsnd — 0.90, 0.95, 0.99 CL allowed regions obtained 
by our analysis for LSND data are also shown, as dark 
grey shaded, light grey shaded, and white areas, respec- 
tively. We find no overlap between the two individual 
95% CL allowed regions; on the other hand, there is over- 
lap between the two 99% CL regions. 

Fig. shows the same (3+1) allowed regions as Fig.|2^ 
but in the (plsnd, Am^) plane, where Plsnd is defined as 
the I'e oscillation probability averaged over the 

LSND L/E distribution: 

Plsnd = {Pi^^t^ v^)) (7) 

where P{i'^ — > Ve) is given by Eq. ^for a = /i, P = e, 
and is a function of all the mass and mixing parameters 
of the oscillation model under consideration. This has 
the obvious disadvantage of being a quantity dependent 
upon the specifics of a certain experiment, as opposed to 
a universal variable such as sin^ 26'^e — ^U^aU^a- How- 
ever, Plsnd bas the advantage of being unambiguously 
defined for any number of neutrino generations, and thus 
is useful in discussing (3+2) models later in this paper. 
As stated previously, the oscillation probability estima- 
tor sin^ 29fj,e = ^Uf^iU'^i cannot be used when more than 
one Am^ value affects the oscillation probability, as is the 
case for (3+2) models. A second advantage of using Plsnd 
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instead of sin^ 261^6 as the oscillation probability estima- 
tor, is that the allowed values for Plsnd inferred from 
the LSND result tend to be almost Am^-independent 
(see grey-shaded areas in Fig. ^p), as expected for an 
almost pure counting experiment such as LSND. The os- 
cillation probability reported by the LSND collaboration 
P is Plsnd = (0.264 ± 0.067 ± 0.045)%, and agrees weh 
with our result of Fig. I^b- 

Fig. |2t shows the values for (xNSBL)mi„ as a function 
of Am^. The number of degrees of freedom is 143. 
As discussed in Appendix ^ the dip in (Xnsbl)™!,, at 
Am^ ~ 0.9 eV^ is due to Bugey data preferring Uei ^ 
values, while the minimum at Am^ 10 — 30 eV^ 
is due to CDHS (mostly) and CCFR84 data, prefer- 
ring Uf^4 ^ values. The value for no oscillations, 
(XNSBL)no™c = 132.2, is the largest value in Fig. [2:; 
this means that the choice of parameters Ud = ?7^4 = 
provides the best-fit to NSBL data, for the Arri^ values 
satisfying the condition (x^sbl)„,„ = (XNSBL)n„ „,c- Note 
that the Am^ ~ 0.9 eV^ Am^ ~ 10 - 30 eV^ dips in 
Xnsbl consistent with Uf.iU^4^ = 0, and therefore with 
sin^ 20^e — Plsnd — 0, but give better fits than the no- 
osciUations hypothesis, J7e4 = C^^4 0. In other words, 
the goodness of fit for the sin^ '2.6 = Plsnd — region 
depends on the Am? value considered. 

Similarly, Fig. [2Ji shows the values for (x?sND)min as a 
function of Am^, used to obtain the LSND allowed re- 
gions drawn in Figs. I^t-, lib- 

We now present a slightly different approach to deter- 
mine the statistical compatibility between the NSBL and 
LSND datasets in (3-1-1) models, which will prove useful 
in comparing the (3-f 1) and (3-1-2) hypotheses. 

In Fig. 13 we show the values for the differences 
^Xnsbl ; ^Xlsnd ^ as well as the corresponding confidence 
levels (5nsbl, <5lsnd, as a function of the LSND oscillation 
probability. The curves are for the set of (3-1-1) models 
with the neutrino mass splitting Am? fixed to the best- 
fit value obtained in a combined NSBL-I-LSND analysis 
(see Section Fill Bll . Am^ = 0.92 eV'^, and mixing matrix 
elements J7^4, UeA treated as free parameters. The value 
for Am? is chosen in this way because it represents to a 
good approximation the value for which one expects the 
best compatibility between the two datasets, as can also 
be seen in Fig. In Fig. we map the (J7e4, U^j.^) 
allowed space into an the (plsnd, Xnsbl): (Plsnd,xLnd) 
spaces. For any given value of Plsnd, the minima for 
the Xnsbl and Xlsnd functions are found in the two {UeA, 
U^i) and one {UeiU^i) free parameters available, respec- 
tively. The process is repeated for several Plsnd values, 
and the collection of these minima for all values of Plsnd 
give the two curves in Fig. The individual confidence 
levels (5nsblj ^lsnd, shown in Fig. are obtained from 
^Xnsbl> ^xLnd in the usual way, by assuming one and 
two free parameters for the LSND and NSBL datasets, 
respectively. 

We now address how to extract areas in parameter 
space of a given combined confidence 6 from two indepen- 
dent experimental constraints, in our case obtained via 



the NSBL and LSND datasets, without assuming statisti- 
cal compatibility a priori. The most straightforward way 
(described, for example, in [l^ll^) is to assign a confi- 
dence level S = (5nsbl(<5lsnd + (1 - <5lsnd)/2) to the over- 
lapping part (if any) between the two separate allowed 
regions in parameter space which are found with the con- 
straint (5nsbl = i^LSND- The extra factor (1 — (5lsnd)/2 is 
due to the fact that the LSND allowed region in the os- 
cillation probability is two-sided, and overlap with the 
NSBL result on the same probability is obtained only for 
downward fluctuations in the LSND result, and not for 
upward ones. 

From Fig. Ob, we find overlapping allowed ranges in 
Plsnd for 1 - (J^sbl = 1 - <5lsnd ^ 2.4%. We Conclude 
that, in (3-1-1) models, the LSND and NSBL datasets are 
incompatible at a combined confidence of (5 ~ 96.4%. In 
our opinion, this value does not support any conclusive 
statements against the statistical compatibility between 
NSBL and LSND data in (3-1-1) models, although it rep- 
resents poor agreement between the two datasets. The 
reader should also refer to Section FV Dl where a different 
method to quantify the compatibility between the NSBL 
and LSND results is discussed. Future short-baseline 
constraints on j/^ — > appearance, as well as on j/^ and 
h'e disappearance, should be able to definitively establish 
whether (3-1-1) models are a viable solution to explain 
the LSND signal. 

B. Combined NSBL+LSND analysis 

The second analysis we perform is a combined 
NSBL-fLSND analysis, with the purpose of obtaining 
the (3-1-1) allowed regions in parameter space, in both 
(sin^ 26'^e, Ato^) and (pLSND,Am^) space. A combined 
analysis of this sort assumes statistically compatible re- 
sults. In Section IIII Al we have shown that the LSND 
and NSBL results are marginally compatible, for (3-1-1) 
models. In the following, we refer to the NSBL-I-LSND 
dataset as the short-baseline (SBL) dataset, and we con- 
struct the x^ function: 

Xsbl — Xnsbl ^ Xlsnd (^) 

where the two contributions Xnsbl and xLnd are now si- 
multaneously minimized with respect to the same set of 
three oscillation parameters Am^, Ue4, U^^. Figs.0^ 
andllt. show the 90%, 95%, and 99% CL 3-dimensional 
allowed regions in (Am^, Ue^^ Utj.4) projected onto 
the (sin^ 20^e, Am^) and (plsnd, Am^) 2-dimensional re- 
gions, respectively, from the combined (3-f 1) analysis of 
SBL data. In this combined analysis, we use the same 
Monte Carlo method described in Section Fill Al We de- 
fine the allowed regions in parameter space by performing 
a global scan. For example, the 95% CL allowed region in 
the three-dimensional space (Ato^, UeA, U^a) is obtained 
by reqmring xLl " (XsBL)mi„ < 7.82, where (XsBL)mi„ is 
the global minimum x^ value. Fig. ^ shows the mini- 
mum XsBL values obtained in the combined fit, as a func- 
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FIG. 3: a)x^ difFerences A^nsblj ^Xlsnd, and b) individual confidence levels ^nsbl, <5lsnd, as a function of the LSND oscillation 
probability Plsnd, for the NSBL and LSND datasets. The curves are for (3+1) models with the neutrino mass splitting Am^ 
fixed to the best-fit value Am'^ = 0.92 eV^ from the combined NSBL+LSND analysis, and variable mixing matrix elements 
Ufj,4, Uei- The solid curves refer to the NSBL dataset, the dotted ones to the LSND dataset. The dashed horizontal lines in 
Fig-Et refer to the 90%, 95%, 99% individual confidence levels, the dotted horizontal line gives the combined confidence level 
S = (5nsbl(5lsnd + (1 — '5lsnd)/2) for which the NSBL and LSND datasets are incompatible. 




FIG. 4: Allowed regions in parameter space from a combined analysis of NSBL and LSND data, in (3+1) models, assuming 
statistical compatibility of the NSBL and LSND datasets. Fig. shows the 90%, 95%, and 99% CL allowed regions in 
(sin^ 29 fie, Am^) space, together with the best-fit point, indicated by the star; b) shows the same allowed regions in (plsnd, Am^) 
space; c) shows the minimum value obtained in the combined analysis as a function of Am^. The number of degrees of 
freedom is 148. 



tion of Ato^. Of course, the Xsbl values shown in Fig.^J: 
for any given /S.m? value are larger than the sum of the 
two contributions x^sblj Xlsndj shown in Figs. |3;,d for 
the same Am^ value, since the latter were separately min- 
imized with respect to the oscillation parameters. Sim- 
ilarly, Figs.|3Ji, andlSja show the projections of the 90%, 



95%, and 99% CL allowed regions in (Am^, Uei, Ufj.4) 
onto the (sin^ 26'ee, Am^) and (sin^ 26^^, Am^) space, re- 
spectively, from the combined (3+1) analysis of SBL 
data. The zenith angle distribution of atmospheric muon 
neutrinos provides a constraint to sin^ ff^j^j that is not in- 
cluded in this SBL analysis; mixings to the right of the 
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FIG. 5: Allowed regions in the parameter spaces relevant for i/e and disappearance from a combined analysis of NSBL and 
LSND data, in (3+1) models, assuming statistical compatibility of the NSBL and LSND datasets. Fig.|SJi shows the 90%, 95%, 
and 99% CL allowed regions in (sin^ 2eee, Am^) space, together with the best-fit point, indicated by the star; b) shows the 
same allowed regions in (sin^ 261 (j/j, Am^ ) space. Mixings to the right of the dashed vertical line in Figl^Jj are excluded at 90% 
CL by atmospheric neutrino results |3^. which are not included in this analysis. 



dashed vertical line in FigElD are excluded at 90% CL 
by atmospheric neutrino results [s^l- The global 
minimum is xiBL=144.9 (148 d.o.f.). This value indi- 
cates an acceptable fit, assuming that t he g oodness-of-flt 
statistic follows the standard p.d.f. j23|; for an alter- 
native goodness-of-fit test, the reader should refer to Sec- 
tion lvBl The individual NSBL and LSND contributions 
to the minimum are x^s3l=137.3 and Xlsnd='''-6j re- 
spectively. This best-fit point corresponds to the mass 
and mixing parameters Am^ — 0.92 eV^, Uei — 0.136, 
[7^4 = 0.205. 



IV. RESULTS FOR (3+2) MODELS 

A. Statistical compatibility between NSBL and 
LSND 

Having introduced the relevant oscillation probability 
formula in Eq. \E\ and the statistical estimator Plsnd 
to compare the NSBL and LSND results in Section 
IIII Al we can now quantitatively address the statistical 
compatibility between the NSBL and LSND datasets 
under the (3+2) hypothesis. 

Ideally, we would like to determine the NSBL upper 
limit for Plsnd , for all possible combinations of the mass 
parameters Am|]^, Am^i. This entails performing a 
scan equivalent to the one described in the (3+1) case 
as a function of Am|j, shown in Fig. [3 In practice, 
the CPU-time requirements to pursue this route were 



prohibitive. 

An easier problem to tackle is to determine the 
statistical compatibility between the NSBL and LSND 
datasets only for the (3+2) models with mass splittings 
Amli, Amli fixed to their best-fit values, as obtained 
by the combined NSBL+LSND analysis that we present 
in Section IIVBI In sections IIII Al and IIII Bl we have 
demonstrated that, at least for (3+1) models, this choice 
is a good approximation for the best possible statistical 
compatibility (see Figs. Inland 01 . 

In Fig. we show the behavior of the values 
Axnsbl ^rid Axlsndj ^rid of the confidence levels (5nsbl 
and (5lsnd, as a function of Plsnd, for the set of (3+2) 
models satisfying the requirements Amli = 0.92 eV^, 
Am^j^ — 22 eV^. By analogy with Fig. O we map the 
four-dimensional space (C/e4, Ufj_4, Ue5, U^^) into the two- 
dimensional spaces (Plsnd,Xnsbl) and (Plsnd, xLbl): 
and we plot the minimum x^ values obtained for any 
given value of Plsnd- The confidence levels shown 
in Fig. EJd are obtained from Fig. considering the 
four free parameters (J7e4, U^a, C/e5, C/^s) in the Xnsbl 
minimization process, and the two free parameters 

{UeAU^^A, Ue5Uf,5) for Xlsnd ■ 

From Fig. we find that, in (3+2) mod- 

els, the NSBL and LSND datasets are incom- 
patible at an individual confidence level of 
<^NSBL = <5lsnd = 1 — 0.215 = 78.5%, and at a combined 
confidence level (5 = 5nsbl('5lsnd + (1-'5lsnd)/2) = 70.0%. 
Fig. Elshould be compared to Fig. 13 obtained for (3+1) 
models. A detailed comparison of the (3+1) and (3+2) 
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FIG. 6: a)x differences Axnsbl and Axlsnd) and b) individual confidence levels iSnsbl and (Jlsnd, as a function of the LSND 
oscillation probability Plsnd, for the NSBL and LSND datasets. The curves are for (3+2) models with the neutrino mass 
splittings Amli and Amli, fixed to the best-fit values Amli = 0.92 eV^, Amli = 22 eV^ from the combined NSBL+LSND 
analysis, and variable mixing matrix elements Uei, UfiA, Ue5, J/^s. The solid curves refer to the NSBL dataset, the dotted ones 
to the LSND dataset. The dashed horizontal lines in Fig.jSJs refer to the 90%, 95%, 99% individual confidence levels, the dotted 
horizontal line gives the combined confidence level S = (Snsbl'Jlsnd for which the NSBL and LSND datasets are incompatible. 



hypotheses is presented in Section Ivl 



B. Combined NSBL+LSND analysis 

We now turn to a combined analysis of the NSBL 
and LSND results in (3+2) models, assuming statisti- 
cal compatibility between the two datasets. The pur- 
pose of this combined analysis is to obtain the allowed 
regions in the mass parameter space (Am|;^, Amg^), re- 
gardless of the simultaneous values for the mixing param- 
eters. Results will be shown for Arngj^ > Am^j^; the case 
Am|]^ > Amgj^ can be obtained by simply interchanging 
Amli with Am^j^. The 95% CL allowed region is de- 

< 



10- 



^SBL 
,2 



(Xsbl)"!*" 



fined as the (Aml^, Am|]^) for which x] 
5.99, where {XsBiJrnin is the absolute minimum for 
all (Am|i,Am2^ ) values. In the minimization proce- 
dure, the mixing matrix elements C/e4, ?7^4, C/e5, C^^ts, 
are treated as free parameters. Fig. [7| shows the 90% 
and 99% CL allowed regions in (Am|j, Amg^) space ob- 
tained in the combined (3+2) analysis. In light of the 
(3+1) analysis shown in previous sections, the result is 
not surprising, pointing to favored masses in the range 
Amli ~ 0.9 eV^ Aml^ -10-40 eV^ at 90% CL. At 
99% CL, the allowed region extends considerably, and 
many other {Amli, Am'^i) combinations appear. The 
best-fit model (xsj3l=135.9, 145 d.o.L) is described by 
the following set of parameters: Ato|]^ — 0.92 eV^, 
Ue4 = 0.121, U^,4 = 0.204, Am^i = 22 eV^ Ue5 = 0.036, 
and Ufj.5 = 0.224. We note here that the best-fit is not 
obtained for fourth and fifth mass eigenstates with de- 
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FIG. 7: Allowed ranges in (Am^i, Amli) space for (3+2) 
models, for the combined NSBL+LSND analysis, assum- 
ing statistical compatibility between the NSBL and LSND 
datasets. The star indicates the best-fit point, the dark and 
light grey-shaded regions indicate the 90 and 99% CL al- 
lowed regions, respectively. Only the Am|i > Am|i region 
is shown; the complementary region Amli > Am^i can be 
obtained by interchanging Am^i with Am§i. 



generate masses, that is for Amli ~ Am^i. The best- fit 



model we found for sub-eV neutrino masses is: 



Am|i 
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0.46 eV^ UeA = 0.090, U^^ = 0.226, Am^^ = 0.89 eV^ 
Ueb = 0.125, Uf^4 = 0.160, corresponding to x^^^^^l^LS 
(145 d.o.f.). 



V. COMPARING THE (3+1) AND (3+2) FITS 
TO SBL DATA 

In this section, we discuss four statistical tests that 
can be used to quantify the better overall agreement of 
SBL data to a (3+2) hypothesis for neutrino oscillations, 
compared to a (3+1) one. 



A. Test 1: NSBL upper limit on Plsnd at a given 
confidence level 5nsbl 

Test 1 uses only NSBL data to establish the (3+1) and 
(3+2) upper bounds on the LSND oscillation probability 
Plsnd- From Figs. |31 and El we obtain at a confidence 
level ^NSBL = 0.90 (0.99): 

(3+1): Plsnd < 0.100% (0.162%) 
(3+2): Plsnd < 0.186% (0.262%) 

Therefore, we find that (3+2) models can enhance the 
LSND probability Plsnd by quite a large factor, com- 
pared to (3+1) models. The increase in Plsnd that 
we obtain is significantly larger than the 25% increase 
found in ,2^], which is based on a specific choice of 
mixing parameters, as opposed to the complete param- 
eter scan performed in this work. The value for the 
— > i?e oscillation probability measured by LSND ^ is 
p^g^j, = (0.264 ± 0.067 ± 0.045)%, where the errors refer 
to the la statistical and systematic errors, respectively. 



B. Test 2: statistical compatibility between the 
NSBL and LSND datasets 

Test 2 uses both the NSBL and LSND datasets, and 
treats them independently to find the combined confi- 
dence level S = (5nsbl(<5lsnd + (1 - <^lsnd)/2) for which 
the datasets are incompatible, both in (3+1) and (3+2) 
models. The combined confidence levels can also be read 
from Figs . 131 and |H1 

(3+1): (5 = 96.4% 
(3+2): 5 = 70.0% 

Therefore, we find that in (3+1) models the two datasets 
are marginally compatible, and the agreement is better 
in (3+2) models. 



C. Test 3: likelihood ratio test 

Test 3 combines the NSBL and LSND datasets into a 
single, joint analysis. The likelihood ratio test [3^ pro- 
vides a standard way to assess whether two hypotheses 
can be distinguished in a statistically significant way. We 
define the maximum likelihood Ci from the minimum 
values (XsBL)mm,i as A = exp(-(xsBL)mm,i/2), where 
the index i = 1, 2 refers to the (3+1) and (3+2) hy- 
potheses, respectively. We can then form the likelihood 
ratio Ai^2 = ^i/^2- If the (3+1) hypothesis were as ad- 
equate as the (3+2) hypothesis in describing SBL data, 
the quantity 

Xl2(3) ^ -21nAi,2 = ixL.Un.i - (xLJ».n,. (9) 

should be distributed as a distribution with three de- 
grees of freedom, where the number of degrees of freedom 
is the difference in the number of mass and mixing pa- 
rameters in the (3+2) and (3+1) hypotheses, 6-3=3. 
In our combined fits, we obtain (see Sections [ill Bl and 

(3+1): (xLl)™,i = 144.9, (148 d.o.f.) 

(3+2): (xLl)™,2 = 135.9, (145 d.o.f.) 

and therefore Xi,2(3) = 9.0. This value is significantly 
larger than 3: the probability for a distribution with 
three degrees of freedom to exceed the value 9.0 is only 
2.9%. In other words, according to the likelihood ratio 
test, the (3+1) hypothesis should be rejected compared 
to the (3+2) one at the 97.1% CL. Therefore, based on 
this test, we conclude from test 3 also that (3+2) models 
fit SBL data significantly better than (3+1) models. 

D. Test 4: compatibility using the "parameter 
goodness-of-fit" 

Test 4 uses both the results of the individual NSBL 
and LSND analyses, as well as the results of the com- 
bined NSBL+LSND analysis. The test is based on the 
"parameter goodness-of-fit" 36] to compare the com- 
patibility of the NSBL and LSND results under the 
(3+1) and (3+2) hypotheses. The test avoids the prob- 
lem that a possible disagreement between the two re- 
sults is diluted by data points which are insensitive 
to the mass and mixing parameters that are common 
to both datasets. The number of parameters com- 
mon to both datasets is Pc = 2 in (3+1) models, and 
Pc = 4 in (3+2) models. One possible choice of com- 
mon parameters is (Aml^, ?7e4J7^4) for (3+1) models, 
(Am|i,J7e4J7^4, Am|i,J7e5;7^5) for (3+2) models. The 
test is based on the statistic xIg = Xpg,nsbl + Xpg,lsnd: 

where XpG.NSBL = (XnSBl)sbL mi„ — (XnSBl)nSBL mln ^ud 
XPG.LSND = (XlSNd)sBL min ~ (XlSNd)lSND min thC (pOS- 

itive) differences for the NSBL and LSND x^ values ob- 
tained by minimizing the entire SBL x^ function, mi- 
nus the x^ values that best fit the individual datasets. 
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Model 


2 

Xpg,nsbl 


2 

Xpg,lsnd 


2 

Xpg 




PG (%) 


(3+1) 


11.8 


4.3 


16.1 


2 


3.2 • 10"^ 


(3+2) 


7.1 


4.4 


11.5 


4 


2.1 



TABLE I: Parameter goodness-of-fit PG, as defined in [s^, 
to test the statistical compatibility between the NSBL and 
LSND datasets under the (3+1) and (3+2) hypotheses. The 
quantities Xpg.nsbl a-nd Xpg.lsnd a^re the NSBL and LSND 
contributions to the test statistic Xpg defined in the text; Pc 
indicates the number of parameters common to both datasets. 



Table gives the values for the parameter goodness- 
of-fit PG as defined in based on the x^g statistic, 
and the number of parameters common to the NSBL and 
LSND datasets, Pc- This test shows a dramatic improve- 
ment in the compatibility between the NSBL and LSND 
results in going from a (3+1) to a (3+2) model, rais- 
ing the compatibility by nearly two orders of magnitude, 
from 0.03% to 2.1%. It will be interesting to investigate if 
(3+3) models can improve the compatibility further. The 
resulting compatibility levels obtained with the param- 
eter goodness-of-fit method are lower than those found 
in Section [VB I this, however, is not surprising, since the 
two statistical tests are quite different. 



VI. ADDITIONAL CONSTRAINTS 

The (3+1) and (3+2) models discussed in this work 
should be confronted with additional experimental con- 
straints, other than the ones discussed in detail in the 
previous sections. We limit ourselves here to list and 
comment on some of these constraints, rather than ad- 
dress them in a quantitative way. Mostly, we will discuss 
the impact that such additional constraints may have on 
the best-fit (3+1) and (3+2) models found in Sections 
IrnlandHVI 

First, nonzero mixing matrix elements C/e4, U^J,'^, Ue5, 
and ?7^5 may cause observable effects in atmospheric neu- 
trino data, in the form of zenith angle-independent sup- 
pressions of the and Ve survival probabilities. Since 
our analysis of SBL data tends to give larger values for 
muon, rather than electron, flavor content in the fourth 
and fifth mass eigenstate, the effect should be larger on 
muon atmospheric neutrinos. For example, the (3+1) 
and (3+2) best-fit models from Sections 1111 Bl and II V Bl 
would give an overall suppression of the flux of 8% 
and 17%, respectively. The size of the effect of Vx 
oscillations at high Am^ is comparable to the current 
accuracy with which the absolute normalization of the 
atmospheric neutrino flux is known |37j . which is ap- 
proximately 20%. A more quantitative analysis using 
the full Sup er-Kamiokande and MACRO spectral infor- 
mation [3JI puts an upper bound of 16% at 90% CL on 
this high Am^ contribution to the atmospheric flux 
suppression (in the notation of Ref. (34] . this suppres- 



sion is parametrized as 2d^(l — d^), where < 0.09 at 
90% CL). Therefore, it is expected that the inclusion of 
atmospheric neutrino data in this analysis would pull the 
best-flt muon flavor components in the fourth and fifth 
mass eigenstates to lower values, but not in a dramatic 
way (see also Fig. ]^). 

Second, models with large masses rui and rns, and 
with nonzero mixing matrix elements Ue4 and Ue5 , should 
be confronted with tritium (3 decay measurements. The 
presence of neutrino masses 7714 and 7715 introduces kinks 
in the differential f3 spectrum; the location in energy of 
the kinks is determined by the neutrino masses, and the 
size of the kinks is determined by the amount of electron 
flavor component in the fourth and fifth mass eigenstates. 
For a spectrometer integrating over the electron energy 
interval 6 near the /3-decay endpoint, the count rate is 

M 

where the quantity R does not depend on the small neu- 
trino masses and mixings, n = 4 or tt, = 5 for (3+1) 
or (3+2) models, respectively, and we have assumed 
d > rrii, i = 1,...,??., and CP invariance. From the 
experimental point of view, tritium /3 decay results are 
generally expressed in terms of a single effective mass 

n,(5) = |(52-m(z.e)^)3/2 (11) 

where ra{ve) is the flt mass parameter. In the limit 6^ ^ 
m^, i = 1, . . . ,n the relation between the true masses and 
mixings to the fltted mass m{Vf.) is independent from the 
integration interval 5: 

n 

The condition ^ mf, i = 1, . . . , n is generally satisfled 
for the neutrino masses considered in this paper, in order 
to ensure sufficient (3 decay count rate statistics in the ex- 
periments. Therefore, to a first approximation, we can 
consider the effect of heavy neutrino masses m^, only 
on the single mass parameter m{ve) fitted by the experi- 
ments. A more general analysis assessing the sensitivity 
of current and future /3 decay experiments to multiple 
fitted neutrino masses, although highly desirable, is be- 
yond the scope of this work; for further details, the reader 
should consult Ref. 38] . The current best measurements 
on m{ve)'^ come from the Troitsk and Mainz experiments 
|39l |. which have very similar m{i'e)^ sensitivities. Both 
found no evidence for a nonzero m{i/f.Y value; the lat- 
est Mainz result is m{vef = -1.6 ± 2.5 ± 2.1 eV^, or 
m{Pe) < 2.2 eV at 95% CL, using S = 70 eV Now, 
assuming a normal hierarchy (mi < 7714 < m^) with 
TOi ~ 0, the (3 decay neutrino mass in Eq. 1121 can be 
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written as m{ve) — U^^Aml^ + U^^Am'^^; the best-fit 
(3+1) and (3+2) models found in this analysis would 
give m(iye)2 = 0.017 eV^ and m(z^e)2 = 0.042 eV^ re- 
spectively, that is m(i^e)^ values well below the current 
experimental sensitivity. The planned tritium /? decay 
experiment KATRIN should be able to improve the sen- 
sitivity to m{ve) by roughly an order of magnitude in the 
forthcoming years, thanks to its better statistics, energy 
resolution, and background rejection Specifically, 
the systematic and statistical (for (5 > 30 eV) uncer- 
tainties on the single fitted mass m{i'e)^ quoted for KA- 
TRIN are 0.018 and 0.016 eV^, respectively ji^, which 
should provide some sensitivity to the preferred (3+1) 
and (3+2) neutrino models with a normal mass hierar- 
chy, mi < TO4 < ms. We now consider mass spectra 
with an inverted hierarchy, defined here as TO4 < mi for 
(3+1) models, and < TO4 < nii for (3+2) models. 
We note that for (3+2) models other hierarchies are also 
possible, but those do not satisfy the implicit assumption 
|Ar7ig]^| = |Am54| + |Am|j^| taken in this analysis. The 
/3 decay neutrino mass in Eq. 1121 can now be written 
as m(i^e)^ — l^'TiliI for inverted (3+1) models, and as 
m(i/e)^ — lAwg]^! for inverted (3+2) models. Clearly, in 
this case the /3 decay constraints depend strongly on the 
values of |Am4i |, | Am|]^ |, and models with > 5 eV^ mass 
splittings are already severely disfavored. 

Third, introducing sterile neutrinos may affect a num- 
ber of cosmological predictions, which are derived from 
various measurements j4l[ . The standard cosmological 
model predicts that sterile neutrinos in the ^ 1 eV mass 
range with a significant mixing with active neutrinos 
would be present in the early Universe with the same 
abundance as the active neutrino species, in disagreement 
with cosmological observations [43 . 143^ . On the other 
hand, several models have been proposed that would rec- 
oncile sterile neutrinos with cosmological observations, 
for example suppressing thermal abundances for sterile 
neutrinos (see, e.g., Ref.j42iJ and references therein). In 
particular, active-sterile oscillations in the early Universe 
may provide a natural mechanism to suppress the relic 
abundances of sterile neutrinos 44] , and scenarios invok- 
ing multiple sterile neutrinos are being investigated |43| . 



VII. CONCLUSIONS 

We have performed a combined analysis of seven short- 
baseline experiments (Bugey, CHOOZ, CCFR84, CDHS, 
KARMEN, LSND, NOMAD) for both the (3+1) and the 
(3+2) neutrino oscillation hypotheses, involving one and 
two sterile neutrinos at high Am^ , respectively. The mo- 
tivation for considering more than one sterile neutrino 
arises from the tension in trying to reconcile, in a CPT- 
conserving, four-neutrino picture, the LSND signal for 
oscillations with the null results obtained by the other 
short-baseline experiments. Multiple {e.g. three) sterile 
neutrinos can also be motivated on theoretical grounds. 

We have described two types of analyses for both the 



(3+1) and (3+2) neutrino oscillation hypotheses. In the 
first analysis, we treat the LSND and the null short- 
baseline (NSBL) datasets separately, and we determine 
the statistical compatibility between the two. In the sec- 
ond analysis, we assume statistical compatibility and we 
combine the two datasets, to obtain the favored regions 
in neutrino mass and mixing parameter space. 

The main results of the analysis are summarized in Sec- 
tion0 where we compare the adequacy of the (3+1) and 
(3+2) hypotheses in describing neutrino short-baseline 
data, by means of four statistical tests. First, we treat 
the LSND oscillation probability as a parameter that can 
be measured with NSBL data alone, and find that the 
NSBL 90% CL upper limit on the LSND oscillation prob- 
ability can be significantly relaxed by going from (3+1) 
to (3+2) models, by about 80%. Second, the combined 
confidence level for which the NSBL and LSND datasets 
are incompatible is determined to be 96.4% and 70.0% in 
the analysis, for the (3+1) and (3+2) hypotheses, respec- 
tively. Third, a likelihood ratio test of the two hypotheses 
is discussed, and shows that the (3+1) hypothesis should 
be rejected compared to the (3+2) one at the 97.1% CL. 
Fourth, the parameter "goodness-of-fit" defined in |3^ 
shows much better agreement between the NSBL and 
LSND results for (3+2) models than for (3+1) models. 

In conclusion, we find that (3+1) models are only 
marginally allowed when considering all of the seven 
short-baseline results, including LSND, in agreement 
with previous analyses H^H^Iloll, and that (3+2) mod- 
els can provide a better description of the data. Only the 
simplest neutrino mass and mixing patterns have been 
fully characterized in the literature so far, and the anal- 
ysis described in this paper may be viewed as a simple 
attempt to explore more generic scenarios, which appear 
both experimentally and theoretically plausible. Given 
the bright potential for precision measurements by neu- 
trino oscillation experiments in the near future, a more 
general phenomenological approach may be needed. 
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APPENDIX A: PHYSICS AND STATISTICAL 
ASSUMPTIONS 

In this section, we briefly describe the physics and sta- 
tistical assumptions used to obtain the approximate char- 
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acterizations of the short-baseline experiments used in 
the analysis. For the analysis of the Bugey, CDHS, and 
KARMEN data, we also refer to the excellent reference 
[3^. which we followed closely. 

The Bugey experiment 18] is sensitive to v,, disappear- 
ance by measuring the charged-current interaction of v^s 
produced by two nuclear reactors at the Bugey nuclear 
power plant. Two liquid scintillator detectors, located at 
different positions, are used. The signature for an an- 
tineutrino interaction is a positron and a delayed light 
pulse produced by the subsequent neutron capture on 
®Li. Data are given for three baselines: 15, 40, and 95 
m between neutrino production and detection. We follow 
the "normalized energy spectra" analysis discussed in the 
Bugey paper 0| . The data are presented as ratios of ob- 
served to predicted (for no oscillations) positron energy 
spectra, between 1 and 6 MeV positron energy. We use 
25,25, and 10 positron energy bins for the 15, 40, 95 m 
baselines, respectively. In the analysis, fits included 
not only the mass and mixing parameters, but also five 
large scale deformations of the positron spectrum due 
to systematic effects. The experimental positron energy 
resolution and the neutrino baseline smearing are taken 
into account; the neutrino cross-section energy depen- 
dence within a positron energy bin is not (the energy bin 
widths are small). 

Similarly, the CHOOZ experiment investigates Pe 
disappearance by observing interactions of 's produced 
by two nuclear reactors ~ 1 km away from the CHOOZ 
detector. The signature for a neutrino interaction is a 
delayed coincidence between the prompt signal and 
the signal due to the neutron capture in the Gd-loaded 
scintillator. We follow "analysis A," as discussed in the 
CHOOZ paper [23|. Data are given as positron yields 
as a function of energy. In this analysis, seven positron 
energy bins, between 0.8 and 6.4 MeV, are considered, 
for which the CHOOZ observations, as well as the pre- 
dictions on the positron yields for the no-oscillation case 
from both reactors, are given in Because of the pres- 
ence of two reactor sources, the analysis comprises 14 
positron yield bins for a given energy/baseline. We use 
the full covariance matrix to take into account the fact 
that the yields corresponding to the same energy bin are 
extracted for both reactors simultaneously, as is done in 
|2ll |. The analysis fits for the systematic uncertainty in 
the absolute normalization constant on the Vf, yield from 
the reactors, in addition to the mass and mixing parame- 
ters. Since we are interested in the Am^ > 0.1 eV^ range 
only, where no energy shape distortions are expected, we 
neglect the systematic uncertainty on the energy-scale 
calibration, and the effect of the positron energy resolu- 
tion. 

The CCFR84 experiment constrains and 
disappearance by measuring the charged-current inter- 
action of muon neutrinos and antineutrinos, produced 
by a Fermilab secondary, sign selected beam yielding 
< El, < 230 GeV neutrinos from tt^ and de- 
cays in the 352 m long decay pipe. We refer here to 



the 1984 CCFR experiment (hence the label CCFR84 
throughout the text), which operated with two similar 
detectors located at different distances from the neutrino 
source, 715 and 1116 m from the mid-point of the decay 
region, respectively. The two sampling calorimeter de- 
tectors consisted of steel plates and scintillation counters. 
Six secondary beam momentum settings were used: five 
for neutrino running, and one for antineutrino running. 
For each secondary beam momentum setting, the data 
are divided into three neutrino energy bins, for a total of 
eighteen energy bins, from Ref . [45l| . Data are presented 
as double ratios: the far to near detector ratio of ob- 
served number of events, divided by the far to near ratio 
of events predicted for no oscillations. As in 19], only 
the mean neutrino energy for a given neutrino energy bin 
is used in the analysis. The systematic and statistical 
uncertainties on the far to near ratio normalization are 
taken into account. The systematic uncertainty is as- 
sumed to be energy-independent and totally correlated 
between any two energy bins. The neutrino pathlength 
smearing, mostly due to the long decay region, is also 
taken into account. 

The CDHS experiment is also sensitive to v^^ dis- 
appearance via the charged-current interaction of I'^i's, 
produced by a 19.2 GeV/c proton beam from the CERN 
Proton Synchrotron. Two detectors are located at 130 
and 835 m from the target. The detectors are sampling 
calorimeters, with iron and scintillator modules inter- 
spersed, to measure the range of a muon produced in a 
neutrino interaction. Fifteen muon range bins are used. 
The data are presented as double ratios: the far to near 
detector ratio of the observed number of events, divided 
by the far to near ratio of the number of events pre- 
dicted for no oscillations. Neutrino energy distributions 
are obtained for a given muon energy (or range) via the 
NUANCE neutrino cross-section generator. As for 
CCFR84, the systematic uncertainty on the far to near 
ratio and the neutrino baseline smearing are taken into 
account. 

The KARMEN experiment investigates the — > 
i>e appearance channel, from P^'s produced in the 7r+- 
/i '''-decay at rest (DAR) chain of the ISIS neutrino 
source. KARMEN measures the charged-current inter- 
action p(Pe,e"'")n, with a liquid scinitillator detector lo- 
cated at an average distance 17.7 m downstream of the 
neutrino source. The signature is a spatially correlated 
delayed coincidence between a prompt positron and a de- 
layed 7 event from a {n, 7) neutron capture reaction. In 
this analysis, only the positron ("prompt") energy dis- 
tribution after all cuts is taken into account, given in 
|22|. The data are binned into nine prompt energy bins, 
between 16 and 50 MeV (all bins are 4 MeV wide, ex- 
cept the highest energy one, ranging from 48 to 50 MeV). 
In predicting the prompt energy distribution for a set of 
mass and mixing oscillation parameters, the given Monte 
Carlo positron energy distribution, and the total num- 
ber of events expected after all cuts for full mixing and 
Am^ — 100 eV^, are used ^3- Energy resolution and 
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FIG. 8: 90% CL upper limits on oscillations derived in this analysis for the following NSBL experiments: a) Bugey, b) CHOOZ, 
c) CCFR84, d) CDHS, e) KARMEN, g) NOMAD. Fig. H shows the LSND 90% CL allowed region obtained with the decay-at- 
rest analysis described in Appendix^ (solid line), superimposed to the published LSND 90% CL allowed region (dashed line). 
Also shown are the (x^)miii values as a function of Am^ obtained by all the experiments considered individually. The number 
of degrees of freedom is 58 in Bugey, 12 in CHOOZ, 16 in CCFR84, 13 in CDHS, 7 in KARMEN, 3 in LSND, 28 in NOMAD. 



baseline smearing effects (due to finite detector size) are 
taken into account. Given the low statistics of the nine 
KARMEN prompt energy bins, we construct the ^ func- 
tion by first defining the Hkelihood ratio 25j : 



/(n;/i(g),b) 
/(n;n,b) 



where Q denotes schematically all mass and mixing pa- 
rameters, n, and b are the data, expected signal, 
and expected background vectors with nine elements, and 
/(n; /i(6'), b) are the probabilities for a Poisson process 



with known background: 

9 



/(n;,(^),b)^n ^'^^'^^"'';^^^''^'^^^ (A2) 



(Al) We define as- 

VVC ucillic A KARMEN 



X 



KARMEN 



-21nA(6') 



(A3) 



The LSND experiment at Los Alamos Q is also 
sensitive to — > Df, appearance, with a neutrino source 
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and detection signature similar to that of KARMEN, 
but with better statistics. The LSND hquid scintillator 
detector is located at an average distance of 30 m from 
the neutrino source. As for KARMEN, in this analysis 
we consider only the positron energy distribution arising 
from a Ve interaction in mineral oil, published as five 
energy bins between 20 and 60 MeV [g. Our analysis 
ignores the information arising from the higher-energy 
neutrinos from pions decaying in flight, which has a 
smaller (but non-negligible) sensitivity to oscillations 
compared to the decay at rest (DAR) sample considered 
here. In our simulation, we take into account the 
expected energy distribution from /i+ decay at rest, 
the neutrino baseline distribution for the 8 m long 
cylindrical detector, the neutrino energy dependence of 
the cross-section for the detection process p(Pe,e+)n 
(including nuclear effects, simulated with the NUANCE 
|4f| neutrino cross-section generator), and the experi- 
mental energy resolution. We use the published numbers 
for the background expectations, the number of i>f. 
events for 100% — > transmutation, and for the 
efficiency of the event selection criteria. We construct 
the LSND function in the same way as we construct 
the one for KARMEN, because of the low statistics of 
the data sample. 

Finally, the NOMAD experiment is sensitive to 
— > Ve oscillations at Am^ > 1 eV^ by looking for 
charged-current muon neutrino and electron neutrino 
interactions in the NOMAD detector |2^. The detector 
consists of a large dipole magnet which houses drift 
chambers to measure the momenta of the charged 
particles produced in neutrino interactions; transition 
radiation modules for lepton identification; an electro- 
magnetic calorimeter to measure the energy of electrons 
and photons; a hadron calorimeter for particle identi- 
fication; and muon chambers for muon identification. 
Neutrinos are produced by impinging 450 GeV protons 
extracted from the CERN SPS accelerator onto a thick 
beryllium target. The secondary particles produced in 
the target are focused into a nearly parallel beam by 
two magnetic lenses, and decay in a 290 m long decay 
tunnel to produce a ~ 10 — 100 GeV neutrino beam with 
about 1% Vf. contamination. Neutrino interactions are 
then observed in the NOMAD detector at an average 
distance of 625 m from the neutrino source. The —>■ v,, 
search is performed by comparing the measured ratio 
Ref_i of the number of Ve to charged-current neutrino 
interactions with the one expected in the absence of 
oscillations. The data are binned into 30 bins, covering 
ten bins in visible energy between 3 and 170 GeV, and 
three radial bins in the neutrino interaction vertex. A 
analysis is performed, using the final NOMAD numbers 
on the observed and predicted electron-to-muon ratio, 
including statistical errors as well as the full error matrix 
describing systematic uncertainties and uncertainty 
correlations over different bins 48]. In predicting the 
effect of — > Ve oscillations under any mass and mixing 
hypothesis, the contribution to i?e^ from oscillations 



with full mixing and Am? = 5, OOP eV^ expected in 
NOMAD after all cuts is used [4g. Energy resolution 
and baseline smearing effects (due to the long decay 
region) are taken into account 23] . 

In Fig. |S1 we show our calculations of the 90% CL 
upper limits on oscillations as a function of Am^ for the 
six NSBL experiments considered here, as well as the 
90% CL allowed region for LSND. The (x^)m.„ values as 
a function of Am^ for all of the experiments are also 
shown. All the solid curves shown are obtained from 



the simplified analysis described here, and compare 


with the published results 


i 


EITi 


m 


m 







The LSND region obtained in our analysis of DAR 
neutrinos is slightly shifted to the right compared to 
the final LSND area, shown in Fig. |Sf as a dashed 
line, reflecting the difference in the two datasets. More 
detailed LSND DAR analyses give results in rough 
agreement with our allowed region 126^ . 

The (x^)min values obtained for the Bugey and CDHS 
experiments as a function of Am? give details that might 
seem surprising, at first. Slightly better fits to the data 
are obtained under a neutrino oscillations hypothesis, as 
opposed to the no oscillations one. Therefore, we add a 
final comment to explain the results of these fits. 

The Bugey fit is driven by the data at the shortest 
baseline, 15 m, where the statistical errors on the 
observed positron spectrum from 9^ interactions are 
the smallest. As explained in Ref . [l^ . systematic 
uncertainties are taken into account by allowing for 
linear deformations, as a function of positron energy, of 
the ratio of observed to predicted positron yields. The 
values of (xBugoy)min as a function of Am^ are explained 
by the fact that, for certain Am^ values, an oscillatory 
fit to the 15 m positron spectrum ratio describes the 
data marginally better than any straight line. Our 
best-fit oscillation hypothesis to Bugey data only is 
Am2 = 0.92 eV^ sin^ 26'ee = 0.05. 

For CDHS, the (XcDHs)mi„ curve in Fig. ISJi has a 
minimum at Am^ ~ 20 — 30 eV^. This minimum 
is due to the fact that the far/near rate ratio, 
corrected for the baseline and detector mass differences 
between the two detectors (as well as other minor 
effects), is measured to be slightly greater than one |20|: 
i?,„„ = 1.044 ± 0.023 ± 0.025. This marginal deviation 
from one causes the fit procedure to prefer more 
disappearance by oscillations in the near than in the 
far detector. Given the average t'^ energy (3.2 GeV) 
and pathlength (130 m) for neutrinos interacting in 
the CDHS near detector, this condition is satisfied in 
the Am^ = 20 — 30 eV^ range. Our best-fit oscillation 
hypothesis to CDHS data only is Am^ = 24 eV^, 
sin^ 20„„ = 0.29. 
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